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Abstract 

In this paper, the familiar problems of free-fall motion and simple harmonic motion (SHM) are 
combined. The novel composite system passes from regular to chaotic behavior for increasing 
values of energy E. This system is a suitable example to introduce undergraduates to some of the 
concepts of dynamical chaos theory. 
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I. INTRODUCTION 



For many centuries now, students of introductory physics have studied Newton's Laws 
of Motion. The two systems done ad nausea are a particle in free-fall and the mass-spring 
system (inclines being a special case of the former). Additionally, it is well know in the chaos 
community, over the past few decades, that coupled oscillators can exhibit chaotic behavior 
i.e. motion that is very sensitive to initial conditions. 1 Examples of this kind include the 
driven damped oscillator, the double pendulum, and the magnetic pendulum. A few of these 
are illustrated in Marion and Thornton's text on classical dynamics. 2 

The purpose of this paper is to answer the question, "What is the behavior of a mass- 
spring system which is subject to free-fall motion above a rigid floor without damping?". 
This system is a variant of a class of problems know as impact oscillators. One of the earliest 
was put forward by Fermi to explain the origin of cosmic rays. 3 His model was an elastic 
ball which bounced off of a vertically oscillating floor. The interested reader is referred to 
the bibliography for descriptions of similar systems.— ^ 

The free-fall spring subject to impact (FSI) constitutes a coupled oscillator problem and 
the resulting motion is chaotic. 



II. THE PHYSICAL SYSTEM: COUPLED OSCILLATORS 

Consider the motion of a point particle released from rest at an initial height y subject 
to the acceleration due to gravity g. If the particle suffers a perfectly elastic collision with 
the floor (as is done in billiard theory) it will return to height y after a time Tf ree -faii given 
by 

T f ree - fall = (1) 

If the system is allowed to continually bounce, then it behaves like an oscillator with 
frequency given by 

^free-fall = Tf, ( 2 ) 

-L free— fall 

Next, consider two equal masses at either end of a perfect spring with an equilibrium 
position of y® and Hooke's constant k. The motion of the difference coordinate i.e. Y = 



2 



2/2 — 2/i will be SHM with frequency 




(3) 



The motion of the two systems can be coupled by releasing the spring system vertically 
and having the lower mass reflect off of a perfectly elastic surface. The system is illustrated 
in figure ^ and represents a coupled oscillator. Examples in nature of such a system would 
be a diatomic molecule of gas which suffers collisions with the floor of its container or a child 
on a pogo stick. 

III. SIMULATION: THE ROUTE FROM PHYSICS TO COMPUTER 

In this section, the steps to get from physical model to a computer simulation are outlined. 
Dimensionless coordinates are also introduced. 

A. Equations of Motion:Model to Theory 

Newton's second law is used to write down the equations of motion. Taking the rigid floor 
to be at y = 0, ni\ = vri2 = m, and the rest length of the spring to be yo the accelerations 
of the respective particles are 



The conventional dot notation has been used to represent the temporal derivatives and 
co>o 2 = k/m . To simplify the above equations EJ dimensionless coordinates are used 




(4a) 
(4b) 
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(5a) 
(5b) 
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(5c) 



J/o^o 2 



The dimensionless quantity (3 contains all of the physical parameters of the problem. 
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The two second order equations can be rewritten into a system of four first order equations 
by noting that velocity is the time derivative of position. In this way, and using system 
the dimensionless equations of motion for the FSI are 

Vi = vi (6a) 
vi = -P + (y 2 - yi -l) (6b) 

y% = v 2 (6c) 

V2 = -(3- (y 2 -yi-1) (6d) 

The above system of equations ignores the interaction of the lower mass yi with the rigid 
floor. The lower mass receives a delta function impulse which switches its velocity v\ — » —V\ 
upon impact. Care must be taken on handling this impact and will be discussed in the 
numerics section. 



B. Energy Bounds 

There are energy bounds that must be taken into consideration for the free-fall spring 
subject to impact (FSI). First, the energy is written down in dimensionless form 

E = ^(V + v 2 2 ) + ^(y 2 — 2/i — l) 2 + y 2 + yi (7) 

The terms are in the order of kinetic, spring potential, and gravitational potential energy. 
The minimum of energy is chosen to be the system at rest with the lower mass on the floor 
at 2/1 = and the upper mass to be at y 2 = 1 — (3 which corresponds to 

E-min = 1 — ^ (8) 

Since negative energy is unphysical in this problem, this relation requires that (3 < 2. 
Moreover, (3 must be chosen so that y 2 > i.e. (3 < 1. By choosing the spring to start off 
at its equilibrium length, the minimum of energy is E m i n = 1. 

The upper bound on energy is found by requiring that the spring not be compressed to 
zero length i.e. y\ — y 2 — 0. This is to ensure that particle two does not suffer a collision 
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with the floor and that the two particles don't pass through each other. Substituting these 
initial conditions into [7| gives 



'max 




(9) 



The energy and parameter requirements for the FSI are 




(10b) 



(10a) 



C. Theory to Numerics 

The system was integrated in time with a fourth order Runge-Kutta method. The time 
step chosen was h = 0.001 which results in a local truncation error of approximately 0(h 5 ) = 
10~ 15 . The energy of the system was monitored and noted to be conserved to eight decimal 
places in double precision arithmetic. For all the simulations (3 was chosen to be (3 = 0.25 
making the energy range 1 < E < 2. 

The subtle part of the simulation is the actual impact of the lower mass yi with the floor 
at y = 0. The time of impact can not be predetermined with a static discrete time step. As 
an example, let yi(U) equal the location of the lower mass right before impact at time ij. At 
the next time step yi(tj + h) , the lower mass has moved through the floor and y\(ti +h) < 
which is unphysical and a consequence of using discrete time steps. The situation can be 
remedied by switching the independent variable from time to y\ and integrating the system 
back to the surface. Applying the chain rule from calculus, the system is transformed to 



dyi 
dyi 



(11a) 



dv\ v\ 
dyi vi 
dy2 _ 2/2 
dyi vi 
dv 2 _ f 2 
dyi vi 



(lib) 



(lid) 



(11c) 
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and integrated with h replaced by h — > — y±. The system , in this way, is moved back 
by the amount of the overshoot. At this point, the velocity of the lower particle is switched 
V\ — > —V\ and the system is integrated forward again with the original time step. A similar 
example can be found in the references as well as a description of Runge-Kutta routines-. 

IV. THE INDICATORS OF CHAOS 

In this section, the system is investigated for the existence of chaos and two popular chaos 
indicators are described. 

A. Sensitivity to Initial Conditions: Lyapunov Exponents 

Chaos has been described as a system which exhibits strong sensitivity to initial condi- 
tions. The famous "butterfly effect" of Edward Lorenz is an example. 7 In figure l2~al two time 
series of yi are plotted on the same graph. Trajectory 1 is at an energy of E — 1.1 and tra- 
jectory two at E — 1.1 + SE where SE = 0.001. The small separation in the two trajectories 
is slightly noticeable and the two paths follow each other closely. Also, the trajectory looks 
periodic. When the energy is increased to E = 1.75 and E = 1.75 + SE something very 
different happens. The two paths in figure 12 .bl start out following each other but around 
t = 10 they start to diverge apart. Also, the motion of neither looks periodic implying that 
higher harmonics have been introduced. The proper way to quantify this phenomena is with 
the Lyapunov characteristic exponent. 

If a system is chaotic, then two trajectories separated slightly in phase space will, on 
average, diverge away from each other with an exponential behavior given by 

d(t) « d e xt (12) 

Where do is the initial phase space separation at t = and A is known as the Lyapunov 
characteristic exponent or the largest Lyapunov exponent. Equation El should not be taken 
too literally but only represents an average behavior of dynamical systems. The calculation 
of A is actually an iterative process. 

To start off, two trajectories initially separated in phase space by an amount do are 
integrate forward for some time r. The system is now separated in phase space by an 
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amount d\. At this point the two trajectories are renormalized so that the separation is 
again do and then integrated forward again for r. Continuing in this manner, a set of N di 
's is used to calculate the Nth approximate of the Lyapunov exponent, given by 



The largest Lyapunov exponent is then given by taking the limit as iV — > oo of equation 
[TBI The interested reader should follow the references for a more detailed description . 2 ! 8 ' 9 ! 10 ' 11 

In the presence of chaos, A tends toward some positive finite value. In the absence of 
chaos, A —>■ 0. In figure EH the results of A for the above two values of energy separated in 
phase space by d = are shown. After 10 6 iterations the results were Xe=i.i ~ 0.004 
and Xe=i.75 ~ 0.187. The results imply that for the smaller energy the motion of the FSI is 
regular and for the larger energy the motion is chaotic. 

B. Poincare maps 

The system under consideration is four dimensional in phase space given by the set 
(yi, 2/2, Vi, Vz). Since energy is a conserved quantity the degrees of freedom can be reduced 
to three. A trajectory through three dimensional phase space of the system would be 
impossible to visualize and too complicated to follow. Information about the system is still 
possible by cutting phase space with a plane and looking at the crossing through this plane 
by trajectories at the same energy. The pattern left behind on the plane is known as a 
Poincare map or a Poincare surface of section.— 

In figure four, four maps for various values of energy are shown. The surface of section 
plane was chosen by plotting points (t>2, 2/2) whenever yi made an impact with the floor. For 
the lowest value of energy, E = 1.1, the motion is regular and periodic. As the energy is 
increased to E = 1.75 the motion becomes more complex and the map becomes dense with 
points. This phase space filling behavior is indicative of chaos. 

V. CONCLUSION 

The free-fall spring subject to impact is a simple model which exhibits rich behavior. 
The reader is invited to carry out other investigations such as higher energies or setting 




(13) 
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mi 7^ m2. The interested reader should see the references for non-technica l 12 i 13 and technical 
introductions to chaos.— 
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The Physical System 
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FIG. 1: Two masses at either end of spring released vertically in a constant gravitational field over 
a perfectly reflecting surface 
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Stable Region: E = 1.1 
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FIG. 4. a: Energy = 1.1 
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FIG. 4.b: Energy = 1.25 
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FIG. 4.c: Energy = 1.5 
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FIG. 4.d: Energy = 1.75 
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